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q^ . Abstract 

O ■ 

c/2 I The adsorption of flexible and highly charged polyelectrolytes onto oppositely charged planar 

ld ' 

2 ' surfaces is investigated by means of Monte Carlo simulations. The effect of image forces stemming 

from the dielectric discontinuity at the substrate interface is considered. The influence, at fixed 

polyelectrolyte volume fraction, of chain length and surface-charge density is also considered. A 

detailed structural study, including monomer and fluid charge distributions, is provided. It is 

demonstrated that image forces can considerably reduce the degree of polyelectrolyte adsorption 

and concomitantly inhibit the charge inversion of the substrate by polyelectrolytes. 
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I. INTRODUCTION 



The adsorption of charged polymers [polyelectrolytes (PEs)] on charged surfaces is an im- 
portant phenomenon in industrial and biological processes. Recent model and well controlled 
experiments |l|,|^ were devoted to characterize PE adsorptions. The understanding of PE ad- 
sorption remains an outstanding problem because of the many different typical interactions 
involved there: strong electrostatic substrate-PE binding, monomer-monomer (PE-PE) re- 
pulsion, chain-entropy, excluded volume, etc. Another complication arises from the dielectric 
discontinuity between the solvent and the substrate generating surface-polarization charges, 
which occurs in experimental situations. 
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On the theoretical side, PE adsorption on 
studied by several authors |3|, y, |5|, la, |7|, ^ 

l23l . |2J| on the level of mean field theories. The case of PE adsorption on heterogeneously 
charged surfaces was recently theoretically addressed by de Vries et al. i21|. A remarkable 
common feature of some of these studies is the charge reversal { overcharQinq) of the substrate 
b. .he adsorbed PE. (see e. g. Refs. QHflflyQ). The PE adsorpfo. onto 
similarly charged substrates were recently investigated by Dobrynin and Rubinstein [19| 
and Cheng and Lai 23|, |2^ . In the latter situation, the PE adsorption is then driven either 
by non- electrostatic short range forces [l9| or attractive image forces |23l . |24| stemming 
from a high-dielectric surface. The problem of repulsive image forces stemming from a low- 
dielectric surface was studied by Borisov et al. |9j and Netz and Joanny [16] on the level of 
the Debye-Hiickel approximation. 

As far as computer simulations are concerned, there exist few Monte Carlo (MC) studies 
about PE adsorption on planar charged substrates 23|, |2y, |27|, |28|, |29|, |30|. The first MC 
study on PE adsorption was that of Beltan et al. |2g| where a lattice-model was employed. 
Yamakov et al. J28| performed extensive MC simulations and found excellent agreement 
with the scaling predictions of Borisov et al. [2|, where different regimes of adsorption 
are identified. Ellis et al. |29'] considered the interesting case of heterogeneously charged 
surfaces (made of positively and negatively charged surface-sites) and demonstrated that a 
PE carrying the same sign of charge as that of the net charge of the substrate can adsorb. 



Cheng et al. 



also investigated the effect of image charges on a high-dielectric constant 



substrate. It is to mention that all these MC simulations l2r 



sibaM 



use the Debye- 



Hiickel approximation. It is only for the case of PE multilayering studied by Messina [30| 
that unscreened long-range electrostatic interactions were properly treated. 

In this article, we investigate multi-chain adsorption in the dilute regime but at fixed PE 
volume fraction in a salt-free environment. Counterions from the substrate and the PEs as 
well as image forces are explicitly taken into account. The influence of chain length (for short 
chains) and substrate-charge density is also considered. Our paper is organized as follows: 
The model and simulation technique are detailed in Sec. m Our results are presented in 
Sec. mil and Sec. IIVI provides concluding remarks. 

II. MODEL AND PARAMETERS 

A. Simulation model 

The setup of the system under consideration is similar to that recently investigated with 
a planar substrate (without image forces) J30l]. Within the framework of the primitive model 
we consider a PE solution near a charged hard wall with an implicit solvent (water at z > 0) 
of relative dielectric permittivity e^ ^ 80. The substrate located at z < is characterized 
by a relative dielectric permittivity e^ which leads to a dielectric jump A^ (when e^ ^ e^) at 
the interface defined as 

A, = ^^^^ > 0. (1) 

The negative bare surface- charge density of the substrate is — o"oe, where e is the (positive) 
elementary charge and ctq > is the number of charges per unit area. Electroneutrality is 
always ensured by the presence of explicit monovalent [Zc = 1) plate's counterions (i.e., 
monovalent cations) of diameter a. PE chains are made up of N^ monovalent positively 
charged monomers [Z^ = 1) of diameter a. Their counterions (monovalent anions) are 
also explicitly taken into account with the same parameters up to to the charge sign as the 
monomers. Hence, all microions are monovalent: Z = Z^ = Z^ = 1 with the same diameter 
size a. 

All these particles making up the system are immersed in a rectangular L x L x t box. 
Periodic boundary conditions are applied in the {x, y) directions, whereas hard walls are 
present at 2; = (location of the charged interface) and z = t (location of an uncharged 
wall) . 



The total energy of interaction of the system can be written as 
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where the first (single) sum stems from the interaction between an ion i (located at z = Zi) 
and the charged plate, and the second (double) sum stems from the pair interaction between 
ions i and j with rij = |rj — rj|. All these contributions to Utot in Eq. Q are described in 
detail below. 

Excluded volume interactions are modeled via a hardcore potential J3l| defined as follows 



Uhs{rij] 
for the microion-microion one, and 



0, for Tij > a 
oo, for Tij < a 



(3) 



jjplate). 



0, for a/2<Zi<T- a/2 
oo, otherwise 



(4) 



for the plate-microion one. For clarity, we recall that a microion stands either for a (charged) 
monomer or a counterion. 

The electrostatic energy of interaction between two microions i and j reads 
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where +(-) applies to microions of the same (opposite) sign, Ib = Pe'^ / iireoer is the Bjerrum 
length corresponding to the distance at which two protonic charges interact with 1//3 = ksT, 
and Ae is given by Eq. (^. The first term in Eq. (0) corresponds to the direct Coulomb 
interaction between real ions, whereas the second term represents the interaction between 
the real ion i and the image of ion j. By symmetry, the latter also describes the interaction 
between the real ion j and the image of ion i yielding an implicit factor 1/2. The electrostatic 
energy of interaction between an ion i and the (uniformly) charged plate reads 
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where, for the first term, +(-) apphes to positively (negatively) charged ions. The second 
term in Eq.® stands for the self-image interaction, i.e., the interaction between the ion i and 
its own image. An appropriate and efficient modified Lekner sum was utilized to compute 
the electrostatic interactions with periodicity in two directions S^- To link our simulation 
parameters to experimental units and room temperature (T = 298K) we choose a = 4.25 
A leading to the Bjerrum length of water Is = 1.68a = 7.14 A. In order to investigate the 
effect of image forces we take a value of e^ = 2 for the dielectric constant of the charged 
substrate (which is a typical value for silica or mica substrates 33|) and e^ = 80 for that of 
the aqueous solvent yielding A^ = |^t| ~ 0.951. The case of identical dielectric constants 
Cp = ^s (Ae = 0) corresponds to the situation where there are no image charges. 

The PE chain connectivity is modeled by employing a standard finite extendible nonlinear 
elastic (FENE) potential for good solvent, which reads 
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with K, = 27kBT/a? and Rq = 1.5a. The excluded volume interaction between chain 
monomers is taken into account via a shifted and truncated Lennard- Jones (LJ) potential 
given by 
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0, for r > l^l^a 

where e = ksT . These parameter values lead to an equilibrium bond length / = 0.98a. 

All the simulation parameters are gathered in Table EJ The set of simulated systems can 
be found in Table |n] The equilibrium properties of our model system were obtained by 



using standard canonical MC simulations following the Metropolis scheme |34l l35 1 . Single- 
particle moves were considered with an acceptance ratio of 30% for the monomers and 50% 
for the counterions. Depending on the parameters, the length of a simulation run ranges 
from 2 X 10^ up to 7 x 10^ MC steps per particle. Typically, about 3 x 10^ to 2.5 x 10^ 
MC steps were required for equilibration, and 1 — 4 x 10^ subsequent MC steps were used 
to perform measurements. 



B. Measured quantities 

We briefly describe the different observables that are going to be measured. In order 
to study the PE adsorption, we compute the monomer density n{z) that is normahzed as 
follows 
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To further characterize the PE adsorption, we also compute the total number of accumulated 
monomers N{z) within a distance z from the planar charged plate that is given by 



N{z) = I n{z')L^dz'. 

'a/2 



(10) 



It is useful to introduce the fraction of adsorbed monomers, N*{z), which is defined as 
follows 

N*{z) = j^. (11) 

Another relevant quantity is the global net fluid charge a(z) which is defined as follows 



a[z] = I [n+{z') - n_{z')]dz', 

la/2 



(12) 



TABLE I: List of key parameters with some fixed values. 



Parameters 



T = 298K 

Ae = or 0. 
Z = 1 
a = 4.25 A 
Ib = 1.68a = 
L = 25a 
T = 75a 

NpE 



951 



7.14 A 



room temperature 

charge number of the substrate 

dielectric discontinuity 

microion valence 

microion diameter 

Bjerrum length 

(x, y)-box length 

z-hox length 

number of PEs 

number of monomers per chain 



TABLE II: Simulated systems' parameters. The number of counterions (cations and anions) 
ensuring the overah electroneutrality of the system is not indicated. 



System 



NpE 



N„ 



<JoL' 



A 
B 
C 
D 
E 
F 
G 
H 



96 
48 
24 
12 
6 

12 
12 
12 



2 

4 

8 

16 

32 

16 

16 

16 



64 
64 
64 
64 
64 
32 
128 
192 



where n+ and «_ stand for the density of all the positive microions (i.e., monomers and 
plate's counterions) and negative microions (i.e., PEs' counterions), respectively. It is useful 
to introduce the reduced surface charge density (y*{z) defined as follows: 



a [z 



aiz 



o-Q 



(13) 



Thereby, cr*{z) corresponds, up to a prefactor aoe, to the net fluid charge per unit area 
(omitting the surface charge density — ctqC of the substrate) within a distance z from the 
charged wall. At the uncharged wall, electroneutrality imposes cr*{z = r — a/2) = 1. By 
simple application of the Gauss' law, [a* {z) — 1] is directly proportional to the mean electric 
field at z. Therefore (J*{z) can measure the screening strength of the substrate by the 
neighboring solute charged species. 



III. RESULTS AND DISCUSSION 

From previous studies pl, UM, 1^ 1^ it is well understood that effects of image charges 
become especially relevant at sufficiently low surface charge density of the interface. It is 
also clear that the self-image interaction {repulsive for A^ > 0, as is presently the case) is 
higher the higher the charge of the ions (polyions) since it scales like Z^. In the present 
situation where we have to deal with PEs, the length of the chain (A^m) is a key parameter 



that can be seen as the valence of a polyion. Hence, we are going to study (i) the influence 
of chain length (Sec. IIIIA|) and (ii) that of surface charge density (Sec. IIIIB|) . For the 
sake of consistency, we fixed the total number of monomers to NpeN^ = 192 meaning 
that the monomer concentration is fixed (see also Table IH^ : The PE volume fraction (p = 

4^ NpMa/2r ^ 2.14 X 10-3 is fixcd. 
3 L^T 

A. Influence of chain length 

In this part, we consider the influence of chain length N^n at fixed surface charge density 
parameter aoL'^ = 64. The latter corresponds experimentally to a moderate [2] surface 
charge density with — o"oe ~ —0.091 C/m^. The chain length is varied from Nm = 2 up to 
Nm = 32 (systems A — E, see Table HTj) . We have ensured that, for the longest chain with 
Nm = 32, size effects are still negligible since the end-to-end distance is about 8.8a which is 
significantly smaller than L = 25a or r = 75a. 

The profiles of the monomer distribution n{z) can be found in Fig. C]and the correspond- 
ing microstructures are sketched in Fig. |21 Let us first comment the more simple case where 
no image charges are present [A^ = - Fig. ^a)]. For (very) short chains (here Nm < 4), 
Fig. ^a) shows that the density profiles exhibit a monotonic behavior even near contact. 
Within this regime of chain length, the monomer density near the charged wall increases 
with increasing A^^- This feature is fully consistent with the idea that stronger lateral cor- 



38l |39| , induce a 



relations, the latter scaling like Z^'"^ for spherical counterions at fixed ctq 
higher polyion adsorption. In other words, at (very) low A^^ conformational entropic effects 
are not dominant and the short-chains systems can be qualitatively understood with the 
picture provided by spherical (or point-like) ions. The scenario becomes qualitatively differ- 
ent at higher chain length [here Nm > 8 - see Fig. ^a)], where n{z) presents a maximum 
near contact which is the signature of a short range repulsion that was also theoretically 
predicted by Borukhov et al. [lO]- This non-trivial feature can be explained in terms of 
entropy: Near the surface of the substrate the number of available PE conformations is 
considerably reduced yielding to an entropic repulsion that can be detected if the driving 
force of PE adsorption (crucially controlled by ao) is not strong enough. This latter state- 
ment will be properly examined and confirmed in Sec. 1111 Bl where the influence of ctq is 
addressed. Nonetheless, the highest value of n{z; Nm) increases with A^ as it should be. All 
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FIG. 1: Profiles of the monomer density n{z) for different chain length N^ with aoL'^ 
(systems A - E). (a) A, = 0. (b) A, = 0.951. 
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these mentioned features can be visualized on the microstructures depicted in Fig. El One 
can summarize those relevant findings, valid for small enough o"o and A^ = 0, as follows: 

• For very short chains the PE adsorption is similar to that occurring with spherical 
electrolytes. 

• PE chains experience a short range repulsion near the substrate due to conformational 
entropic effects. 

We now turn to the more complicated situation where image forces are present [A^ = 
0.951 - Fig. mb)]. An immediate remark that can be drawn from a comparison with the 
Ag = case is that the PE adsorption is much weaker due to the repulsive image-polyion 
interactions. At all Nm-, n{z) presents a maximum dX z = z* that is gradually shifted to 
larger z with increasing A^^- In other words, the thickness of the adsorbed PE layer as 
determined by z* increases with Nm- This phenomenon is of course due to the fact that 
the image-polyion repulsion increases with A^^, similarly to what happens with multivalent 
(point-like or spherical) counterions (3fil . l87| . On the other hand, interestingly, the monomer 
density at contact decreases with increasing A^- This is the result of a combined effect of (i) 
conformational entropy as explained above and (ii) the A^-induced image-polyion repulsion. 
All those features are well illustrated on the microstructures of Fig. |21 
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FIG. 2: Typical equilibrium microstructures of systems A — E. The little counterions are omitted 
for clarity. 
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FIG. 3: Profiles of tlie fraction of adsorbed monomers N*{z) for different chain length Nm (as 
indicated by its numerical value) with (JqL'^ = 64 (systems A — E). (a) A^ = 0. (b) A^ = 0.951. 

To gain further insight into the properties of PE adsorption, we have plotted the fraction 
of adsorbed monomers N*{z) [Eq. (fTT|) ] in Fig. IHl At A^ = [Fig. Efa)], it is observed in the 
immediate vicinity of the wall (roughly for z < 1.5a) that N*{z; Nm) increases monotonically 
with Nm as expected. However, further away from the wall, a non-trivial effect is found 
where N*{z; Nm) surprisingly exhibits a non-monotonic behavior with respect to N^,- More 
explicitly, in the regime of large A^^ we have N*{z;Nm = 32) that is clearly smaller than 
N*{z] Nm = 16) and even smaller than N*{z] Nm = 8) when one is sufficiently far from the 
wall. This remarkable phenomenon is going to be explained later by advocating the role of 
overcharging. Upon switching the image forces on [A^ = 0.951 - Fig. Efb)], N*{z] Nm) shows 
a qualitatively different behavior than that found at A^ = 0, in accordance with our study 
concerning n{z). More precisely: (i) very close to the wall, N*{z;Nm) decreases with A^^ 
whilst (ii) sufficiently far away from the wall N*{z]Nm) increases with A^^. This behavior 
is fully consistent with our mechanisms previously discussed for n{z). Below, we are going 
to show that the reduced net fluid charge (t*{z) is a key observable to account for those 
reported properties of N*{z; Nm)- 

A deeper understanding of the physical mechanisms involved in PE adsorption can be 
gained by considering the net fluid charge parameter cr*{z) [Eq. (fT^] that describes the 
screening of the charged interface. The proffies of (J*{z) for different A^ can be found in 



11 



A =0 



A =0.951 






N 


= 2 





N 


= 4 





N 


= 8 




A',, 


= 16 





W, 


= 32 



4 6 

z/a 



10 




FIG. 4: Profiles of the reduced net fluid charge cr*{z) for different chain length Nm with cioL^ = 64 
(systems A - E). (a) A^ = 0. (b) A^ = 0.951. 

Fig. El At Ae = [see Fig. Efa)], it is shown that for long enough chains (here Nm > 4) 
the substrate gets locally overcharged as signaled by (t*{z) > 1. Physically, this means that 
the global local charge of the adsorbed monomers J4i| is larger in absolute value than that 
of the plate's charge. In other words, the plate is overscreened by the adsorbed PE chains. 
Fig. ^a) indicates that the degree of overcharging increases with N^ as expected from 
the behavior of multivalent counterions, and seems to saturate at high A^^- This enhanced 
Am-overcharging leads to a sufficiently strong effective repulsion between the substrate and 
the PEs in the solution, which in turn prevents from further adsorption. It is precisely 
this mechanism that explains the apparent anomaly found in Fig. Efa) where, sufficiently 
away from the surface, it was reported a significantly lower monomer fraction N*{z; Nm) at 
Nm = 32 than at A^ = 16 or A^^ = 8. This spectacular effect is well illustrated in Fig. 
El (with Nm = 32) where, above the (strongly bound) adsorbed PEs, there is a depletion 
region leading to a plateau in N*{z; Nm = 32). 

Upon inducing polarization charges [A^ = 0.951 - Fig. jlfb)] overscreening is canceled. 
This, in turn, accounts for the absence of plateau in N*{z; Nm) at A^ = 0.951. That striking 
disappearance of overcharging can be rationalized by establishing again an analogy with 
multivalent spherical ions: 

• For the sake of simplicity, let us assume that the PE can be electrostatically envisioned 
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as a spherical polyion of valence Nm with a radius corresponding roughly to the radius 
of gyration of the chain. Thereby, the image-polyion repulsive interactions [including 
the self-image repulsion as well as the lateral image-ion correlations as given by the 
second term of Eq. Q] scale like N'^ whereas the attractive driving force of polyion 
adsorption due to Wigner crystal ordering scales like NJ, |37|. The latter driving 
force corresponds to the highest possible attraction between the substrate and the 
polyion and is therefore a good candidate for the present discussion. Consequently, at 
large enough N^ image forces will be dominant and inhibit overcharging. 



B. Influence of substrate surface-charge density 

To complete our investigation, we would like to address the influence of the substrate 
charge density on the PE adsorption in presence of image forces. In this respect, we consider 
(at fixed Nm = 16) three additional values of the charge density: (TqL'^ = 32, 128, 192 
corresponding to the systems F,G,H, respectively (see Table IIT|l . 

The plots of the monomer density n{z) at various values of ctoL"^ can be found in Fig. 
El Microstructures of systems F and H are presented in Fig. IHl At A^ = [see Fig. 
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FIG. 5: Profiles of the monomer density n{z) for different parameters of surface charge density 
(7oL^ with Nm = 16 (systems D,F — H). The case aoL^ = 64 (system D) from Fig. ^is reported 
here again for easier comparison, (a) A^ = 0. (b) A^ = 0.951. 
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FIG. 6: Typical equilibrium microstructures of systems F and H. The little counterions are 
omitted for clarity. 

Efa)], the monomer density at contact increases with ctq as it should be. Interestingly, the 
local maximum in n{z) [present at small ctq (here o-qL"^ ^ 64)] vanishes at large o"o [see Fig. 
El^a)]. This feature is the result of a cro-enhanced driving force of adsorption that overcomes 
entropic effects at large enough cxo. The strong adsorption at aoL"^ = 192 leads to a flat PE 
layer as well illustrated in Fig. El 

By polarizing the substrate surface (A^ = 0.951), it can be seen from Fig. Efb) and the 
snapshot of Fig. IHl that there is a strong monomer depletion near contact at aoL"^ = 32. 
This feature is due to the combined effects of (i) conformational entropy, (ii) image-monomer 
repulsion and (iii) a lower electrostatic wall-monomer attraction. Upon increasing o"o the 
monomer density near contact becomes larger, and concomitantly, the maximum in n{z) is 
systematically shifted to smaller z. It is to say that the thickness of the adsorbed PE layer 
decreases with ctq- 

The profiles of N*{z) are provided in Fig. [7|from which further characterization of PE 
adsorption can be obtained. At A^ = 0, Fig. [7|^a) indicates that N*{z; o"o) increases with ctq 
but saturates at high ctq. This latter saturation effect should only be relevant for a regime of 
charge where r] = ^^^^2*" is about unity. Indeed, in a typical experimental situation at finite 
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FIG. 7: Profiles of the fraction of adsorbed monomers N*{z) for different parameters of sm'face 
charge density a^L^ (as indicated by its numerical value) with Nm = 16 (systems D, F — H). The 
case aoL^ = 64 (system D) from Fig. El is reported here again for easier comparison, (a) A^ = 0. 
(b) A, = 0.951. 

monomer concentration (even in the dilute regime) we have r] 3^ 1 so that overcharging is 
always possible at large (Tq and thereby N*{z;ao) should always significantly increase with 
o"o as long as packing effects (as generated by the excluded volume of the monomers) are 
not vivid. In parallel, the plateau reported at <7oL^ = 128 and ctq-L^ = 192 in Fig. U^a), 
is the signature of a monomer depletion above the adsorbed PE layer (see also Fig. Ej) due 
to a strong screening of the surface charge by the latter. At A^ = 0.951, Fig. [7|^b) shows 
that N*{z) is considerably smaller than at A^ = even for high ctq, in accordance with the 
behavior of n{z) from Fig. El The A^-induced desorption is especially strong at cqL^ = 32 
where the image-monomer repulsion clearly counterbalances the electrostatic wall-monomer 
attraction. More quantitatively, at z = 3a (a 2;-distance corresponding roughly to the radius 
of gyration of the chain with Nm = 16) about 30% [i.e., N*{z) = 0.3] of the monomers are 
adsorbed with A^ = against only 10% with A^ = 0.951 [see Fig. [Z(b)]. 
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IV. CONCLUDING REMARKS 

We first would like to make some final remarks about the presented results. As far as 
the charge surface distribution on the substrate's surface is concerned, we have assumed a 
smeared- out one in contrast to real experimental situation where it is discrete. Previous 
numerical studies 4l|, |42, |43| have shown that the counterion distribution at inhomoge- 
neously charged substrates may deviate from that obtained at smeared-out ones at strong 
Coulomb coupling (i.e., multivalent counterions and/or high Bjerrum length) or strong sub- 
strate charge modulations. Nonetheless, at standard Bjerrum length (i.e., /^ = 7.1 A for 
water at room temperature as is presently the case) and with discrete monovalent ions 
generating the substrate's surface charge, it has been demonstrated that the counterion dis- 



tribution is marginally modified |4l| even for trivalent counterions. Hence, we think that 
our results will not qualitatively differ from the more realistic situation of non-smeared-out 
substrate charges consisting of discrete monovalent ions. 

An other approximation in our model is the location of the dielectric discontinuity. More 
precisely, it was implicitly assumed that the latter coincides with the charged interface 
(considered here as a hard wall). In fact, experimentally, it is not clear where the dielectric 
discontinuit y is located and the transition is rather gradual and spreads out over several 



Angstroms 4J|, so that in a continuum description the dielectric discontinuity might be 
located somewhat below the hard interface. In this respect our model tends to slightly 
overestimate the effect of image forces and namely, with A^ > 0, the self-image repulsion. 
Furthermore, in the presence of short-range attractive interactions between the substrate 
and the PEs (for instance stemming from some specific chemical properties of the chains 
and the substrate, i.e. chemisorption), the effect of image charges might also be reduced 
[l9| . This means that the substrate-charge undercompensation by PEs induced by repulsive 
image forces as reported in Fig. IH^b) is dependent on the relative strength of that short- 
range attractive interaction, which is not taken into account in our model. Nevertheless, 
we are confident that our results provide a reliable fingerprint for the understanding of the 
effect of image forces on PE adsorption in a salt-free environment. 

It is not a straightforward task to access experimentally to these effects stemming from 
image forces. One major difficulty arises from the fact that by changing the dielectric 
constant of the solvent, ei, one changes the degree of ionization of the PEs. However, there 
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is the experimental possibility to tune A^ by using organic solvents (i.e., with a low ei 
but still polar) with a mixture of large colloidal particles [e.g., latex particles with weak 
curvature and (low) dielectric constant €2 such that €2 < Ci] and PEs. In this experimental 
context, one should be able to verify the trends of our current findings. 

To conclude, we have performed MC simulations to address the effect of image forces on 
PE adsorption at oppositely charged planar substrates. The influence of chain length and 
surface-charge density was also considered. We have considered a finite monomer concentra- 
tion in the dilute regime for relatively short chains. Our main findings can be summarized 
as follows: 

• For very short chains (here Nm < 4) and with no image forces (i.e., A^ = 0), the PE 
adsorption is similar to that occurring with little (spherical) multivalent counterions. 
For longer chains (here Nm > 8), the PEs experience (even at A^ = 0) a short 
range repulsion near the substrate due to chain entropy effects. This latter feature is 
especially relevant at low substrate charge ctq. 

• At fixed ctq and in the presence of repulsive image forces (here A^ = 0.951), it was 
demonstrated that the monomer depletion in the vicinity of the substrate as well as 
the thickness of the PE layer grow with chain length Nm- Concomitantly, the charge 
reversal of the substrate by the adsorbed PEs vanishes. 

• Upon varying gq at fixed N^,, it was shown at A^ = that the net substrate-PE 
force becomes purely attractive at sufficiently high o"o, where chain-entropy effects are 
overcompensated. When image forces are present, the PE depletion near the substrate 
as well as the thickness of the adsorbed PE layer decrease with ctq. 



A future work will address the adsorption of stiff rod-like PEs. This 



very interesting 



22| . Numerical 



situation was recently theoretically examined by Cheng and de la Cruz 

simulation data would be of great help to further characterize the transversal and in-plane 

structures as well as to elucidate the influence of image forces. 
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